PRECAUTIONARY MEASURES FOR CREDIT RISK MANAGEMENT IN JUMP MODELS 



MASAHIKO EGAMI AND KAZUTOSHI YAMAZAKI 

Abstract. Sustaining efficiency and stability by properly controlling the equity to asset ratio is one of 
the most important and difficult challenges in bank management. Due to unexpected and abrupt decline 
of asset values, a bank must closely monitor its net worth as well as market conditions, and one of its 
important concerns is when to raise more capital so as not to violate capital adequacy requirements. In 
this paper, we model the tradeoff between avoiding costs of delay and premature capital raising, and solve 
the corresponding optimal stopping problem. In order to model defaults in a bank's loan/credit business 
portfolios, we represent its net worth by Levy processes, and solve explicitly for the double exponential 
jump diffusion process and for a general spectrally negative Levy process. 
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1. Introduction 

As an aftermath of the recent devastating financial crisis, more sophisticated risk management practices are now 
being required under the Basel II accord. In order to satisfy the capital adequacy requirements, a bank needs to 
closely monitor how much of its asset values has been damaged; it needs to examine whether it maintains sufficient 
equity values or needs to start enhancing its equity to asset ratio by raising more capital and/or selling its assets. 
Due to unexpected sharp declines in asset values as experienced in the fall of 2008, optimally determining when to 
undertake the action is an important and difficult problem. In this paper, we give a new framework for this problem 
and obtain its solutions explicitly. 

We propose an alarm system that determines when a bank needs to start enhancing its own capital ratio. We use 
Levy processes with jumps in order to model defaults in its loan/credit assets and sharp declines in their values 
under unstable market conditions. Because of their negative jumps and the necessity to allow time for completing 
its capital reinforcement plans, early practical action is needed to reduce the risk of violating the capital adequacy 
requirements. On the other hand, there is also a cost of premature undertaking. If the action is taken too quickly, it 
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may run a risk of incurring a large amount of opportunity costs including burgeoning administrative and monitoring 
expenses. In other words, we need to solve this tradeoff in order to implement this alarm system. 

In this paper, we properly quantify the costs of delay and premature undertaking and set a well-defined objective 
function that models this tradeoff. Our problem is to obtain a stopping time that minimizes the objective function. 
We expect that this precautionary measure gives a new framework in risk management. 

1.1. Problem. Let (f2, J 7 , F) be a complete probability space on which a Levy process X = {X t ;t > 0} is 
defined. We represent, by X, a bank's net worth or equity capital allocated to its loan/credit business and model 
the defaults in its credit portfolio in terms of the negative jumps. For example, for a given standard Brownian 
motion B = {B t ; t > 0} and a jump process J = { J t ; t > 0} independent of B, if it admits a decomposition 

X t = x + + aB t + Jt, < t < oo and Xq = x 

for some /i£l and a > 0, then Jt models the defaults as well as rapid increase in capital whereas the non-jump 
terms /xi and oB t represent, respectively, the growth of the capital (through the cash flows from its credit portfolio) 
and its fluctuations caused by non-default events (e.g., change in interest rates). 

Since X is spatially homogeneous, we may assume, without loss of generality, that the first time X reaches or 
goes below zero signifies the event that the net capital requirement is violated. We call this the violation event and 
denote it by 

9 := M{t > : X t < 0} 

where we assume inf = oo. Let F = (Ft)t>o be the filtration generated by X. Then 9 is an F-stopping time 
taking values on [0, oo]. We denote by S the set of all stopping times smaller than or equal to the violation event; 
namely, 

S := {stopping time r : r < 9 a.s.} . 

We only need to consider stopping times in S because the violation event is observable and the game is over 
once it happens. By taking advantage of this, we see that the problem can be reduced to a well-defined optimal 
stopping problem; see Section 2. Our goal is to obtain among S the alarm time that minimizes the two costs we 
describe below. 

The first cost we consider is the risk that the alarm will be triggered at or after the violation event: 



l {r>e,e<oo} 



T£S. 



Here q G [0, oo) is a discount rate and F x is the probability measure and E x ' is the expectation under which the 
process starts at Xq = x. We call this the violation risk. In particular, when q = 0, it can be reduced under a 
suitable condition to the probability of the event {r > 9}; see Section 2. 
The second cost relates to premature undertaking measured by 



H^ h \r) := E x 



l {r<oo} J e~ qt h(Xt)dt 



, t£S. 



We shall call this the regret, and here we assume h : (0, oo) — > [0, oo) to be continuous, non-decreasing and 



(1.1) E a 



e- qt h(X t )dt 



o 



< oo, x > 0. 
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The monotonicity assumption reflects the fact that, if a bank has a higher capital value X, then it naturally has 
better access to high quality assets and hence the opportunity cost h(-) becomes higher accordingly. In particular, 
when h = 1 (i.e., h{x) = 1 for every x > 0), we have 

(1.2) H^(t) = E x - E x t and H^ 1 \t) = - (e*^ - ^] - E^e" 9 "]) , q > 

where the former is well-defined by (1.1). 

Now, using some fixed weight 7 > 0, we consider a linear combination of these two costs described above: 

(1.3) U^ h \r n ) := 4*>(r) + 7 i#' h) (r), r G S. 

We solve the problem of minimizing (1.3) for the double exponential jump diffusion process and a general spec- 
trally negative Levy process. The objective function is finite thanks to the integrability assumption (1.1), and hence 
the problem is well-defined. 

The form of this objective function in (1.3) has an origin from the Bayes risk in mathematical statistics. In the 
Bayesian formulation of change-point detection, the Bayes risk is defined as a linear combination of the expected 
detection delay and false alarm probability. In sequential hypothesis testing, it is a linear combination of the 
expected sample size and misdiagnosis probability. The optimal solutions in these problems are those stopping 
times that minimize the corresponding Bayes risks. Namely, the tradeoff between promptness and accuracy is 
modeled in terms of the Bayes risk. Similarly, in our problem, we model the tradeoff between the violation risk 
and regret by their linear combination U. 

We first consider the double exponential jump diffusion process, a Levy process consisting of a Brownian 
motion and a compound Poisson process with positive and negative exponentially-distributed jumps. We consider 
this classical model as an excellent starting point mainly due to a number of existing analytical properties and 
the fact that the results can potentially be extended to the hyper-exponential jump diffusion model (HEM) and 
more generally to the phase-type Levy model. Due to the memoryless property of the exponential distribution, the 
distributions of the first passage times and overshoots by this process can be obtained explicitly (Kou and Wang 
[20]). It is this property that leads to analytical solutions in various problems that would not be possible for other 
jump processes. Kou and Wang [21] used this process as an underlying asset and obtained a closed form solution 
to the perpetual American option and the Laplace transforms of lookback and barrier options. Sepp [35] derived 
explicit pricing formulas for double-barrier and double-touch options with time-dependent rebates. See also Lipton 
and Sepp [28] for applications of its multi-dimensional version in credit risk. Some of the results for the double- 
exponential jump diffusion process have been extended to the HEM and phase-type models, for example, by Cai 
et al. [7, 8] and Asmussen et al. [1]. 

We then consider a general spectrally negative Levy process, or a Levy process with only negative jumps. 
Because we are interested in defaults, the restriction to negative jumps does not lose much reality in modeling. 
We also see that positive jumps do not have much influence on the solutions. We shall utilize the scale function to 
simplify the problem and obtain analytical solutions. In order to identify the candidate optimal strategy, we shall 
apply continuous and smooth fit for the cases X has paths of bounded and unbounded variation, respectively. The 
scale function is an important tool in most spectrally negative Levy models and can be calculated via algorithms 
such as Surya [40] and Egami and Yamazaki [13]. 
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1.2. Literature review. Our model is original, and, to the best of our knowledge, the objective function defined 
in (1.3) cannot be found elsewhere. It is, however, relevant to the problem, arising in the optimal capital structure 
framework, of determining the endogenous bankruptcy levels. The original diffusion model was first proposed by 
Leland [26] and Leland and Toft [27], and it was extended, via the Wiener-Hopf factorization, to the model with 
jumps by Hilberink and Rogers [17]. Kyprianou and Surya [24] studied the case with a general spectrally negative 
Levy process. In their problems, the continuous and smooth fit principle is a main tool in obtaining the optimal 
bankruptcy levels. Chen and Kou [10] and Dao and Jeanblanc [12], in particular, focus on the double exponential 
jump diffusion case. 

In the insurance literature, as exemplified by the Cramer-Lundberg model, the compound Poisson process is 
commonly used to model the surplus of an insurance firm. Recently, more general forms of jump processes are 
also used (e.g., Huzak et al. [18] and Jang [19]). For generalizations to the spectrally negative Levy model, see 
Avram et al. [2], Kyprianou and Palmowski [23], and Loeffen [29]. The literature also includes computations of 
ruin probabilities and extensions to jumps with heavy-tailed distributions; see Embrechts et al. [14] and references 
therein. See also Schmidli [34] for a survey on stochastic control problems in insurance. 

Mathematical statistics problems as exemplified by sequential testing and change-point detection have a long 
history. It dates back to 1948 when Wald and Wolfowitz [41, 42] used the Bayesian approach and proved the 
optimality of the sequential probability ratio test (SPRT). There are essentially two problems, the Bayesian and 
the variational (or the fixed-error) problems; the former minimizes the Bayes risk while the latter minimizes the 
expected detection delay (or the sample size) subject to a constraint that the error probability is smaller than some 
given threshold. For comprehensive surveys and references, we refer the reader to Peskir and Shiryaev [31] and 
Shiryaev [36]. Our problem was originally motivated by the Bayesian problem. However, it is also possible to 
consider its variational version where the regret needs to be minimized on constraint that the violation risk is 
bounded by some threshold. 

Optimal stopping problems involving jumps (including the discrete-time model) are, in general, analytically 
intractable owing to the difficulty in obtaining the overshoot distribution. This is true in our problem and in the 
literatures introduced above. For example, in sequential testing and change-point detection, explicit solutions can 
be realized only in the Wiener case. For this reason, recent research focuses on obtaining asymptotically optimal 
solutions by utilizing renewal theory; see, for example, Baron and Tartakovsky [3], Baum and Veeravalli [4], Lai 
[25] and Yamazaki [44]. Although we do not address in this paper, asymptotically optimal solutions to our problem 
may be pursued for a more general class of Levy processes via renewal theory. We refer the reader to Gut [16] 
for the overshoot distribution of random walks and Siegmund [37] and Woodroofe [43] for more general cases in 
nonlinear renewal theory. 



1.3. Outline. The rest of the paper is structured as in the following. We first give an optimal stopping model for 
a general Levy process in the next section. Section 3 focuses on the double exponential jump diffusion process 
and solve for the case when h = 1. Section 4 considers the case when the process is a general spectrally negative 
Levy process; we obtain the solution explicitly in terms of the scale function for a general h. We conclude with 
numerical results in Section 5. Long proofs are deferred to the appendix. 
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2. Mathematical Model 

In this section, we first reduce our problem to an optimal stopping problem and illustrate the continuous and 
smooth fit approach to solve it. 

2.1. Reduction to optimal stopping. Fix r £ S, q > and x > 0. The violation risk is 



'1 



{r>e,e<oo} 



E 2 



'1 



{ T =e,e<oo} 



E x [e- qT l {T=e , T<oo} ] 



where the third equality follows because r < 9 a.s. by definition. Moreover, we have 



because 



and 



1{t=0,t<oo} — 1{X t <0,t<oo} a - s - 



{r = 9, t < oo} = {X T < 0, r = 9, r < oo} 



{X T < 0, r < oo} = {X T < 0, r < 6, t < oo} = {X T < 0, r = 9, r < oo} 
where the first equality holds because tG5 and the second equality holds by the definition of 9. Hence, we have 
(2.1) R^(r)=E x [e-< 1T l {XT < Q , T<oo} }, reS. 

For the regret, by the strong Markov property of X at time r, we have 



(2.2) 



H<*>% 



E 2 



l{r<oo} / e- qt h(X t )dt 



E x 



Q^>(X T )l {T<oo} 



where {Q^ q ' h) (X t );t > 0} is an F-adapted Markov process such that 



Q {q ' h \x) :=E 2 



e- qt h{X t )dt 



x € 



Therefore, by (2.1)-(2.2), if we let 

(2.3) G{x) := 1 {1 .< 0} + lQ {q ' h) (x)l {x>0} , i£l 

denote the cost of stopping, we can rewrite the objective function (1.3) as 

C/i^)(T, 7 ) = E x [e- qT G(X T )l {T<oo} ] , reS. 

Our problem is to obtain 

inf E* [e- qT G(X T )l {T<oo} ] 

and an optimal stopping time r* G S that attains it if such a stopping time exists. It is easy to see that G(x) is 
non-decreasing on (0, oo) because h(x) is. If C7(0+) > 1, then clearly 9 is optimal. Therefore, we ignore the 
trivial case and assume throughout this paper that 



(2.4) 



G(0+) < 1. 



As we will see later, when X has paths of unbounded variation, G(0+) = and the assumption above is automat- 
ically satisfied. 
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The problem can be naturally extended to the undiscounted case with q = 0. The integrability assumption (1.1) 
implies K x 9 < oo (without this assumption G(x) = oo and the problem becomes trivial). This also implies 6 < oo 
a.s. and the violation risk reduces to the probability 

R^(t) = ¥ x {t > 9} = F x {t = 6} = F x {X T < 0} , r G S. 

We shall study the case q = for the double exponential jump diffusion process in Section 3. 

2.2. Obtaining optimal strategy via continuous and smooth fit. Similarly to obtaining the optimal bankruptcy 
levels in Leland [26], Leland and Toft [27], Hilberink and Rogers [17] and Kyprianou and Surya [24], the con- 
tinuous and smooth fit principle will be a useful tool in our problem. Focusing on the set of threshold strategies 
defined by the first time the process reaches or goes below some fixed threshold, say A, 

t a := inf{t > : X t < A}, A>0, 

we choose the optimal threshold level that satisfies the continuous or smooth fit condition and then verify the 
optimality of the corresponding strategy. 

Let the expected value corresponding to the threshold strategy t a for fixed A > be 

Mx)--=Ui q ' h) (T A , 7 ), xGR, 

and the difference between the continuation and stopping values be 

(■ta 



(2.5) S A (x) :=4> A (x)-G{x) = R^(t a )- 7 E x I e-*h(X t )dt 

Jo 

We then have 

G(x) + 6 A (x), x > A, 



< A < x. 



(2.6) 4> A {x) 



G{x), 0<x<A, 
1, x<0. 



The continuous and smooth fit conditions are S A (A+) = and S' A (A+) = 0, respectively. For a comprehensive 
account of continuous and smooth fit principle, see Peskir and Shiryaev [31, 32, 33]. 

2.3. Extension to the geometric model. It should be noted that a version of this problem with an exponential 
Levy process Y = {Y t = exp(X t ); t > 0} and a slightly modified violation time 

6 = inf {t > : Y t < a} 

for some a > can be modeled in the same framework. Indeed, defining a shifted Levy process 

X:= \x t = Xi-loga;t>0}, 

we have 

9 = inf jt > : X t < o} . 

Moreover, the regret function can be expressed in terms of X by replacing h(x) with h{x) = h(exp(x + log a)) for 
every x > 0. The continuity and non-decreasing properties remain valid because of the property of the exponential 
function. 
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3. Double Exponential Jump Diffusion 

In this section, we consider the double exponential jump diffusion model that features exponential-type jumps in 
both positive and negative directions. We first summarize the results from Kou and Wang [20] and obtain explicit 
representations of our violation risk and regret. We then find analytically the optimal strategy both when q > 
and when q = 0. We assume throughout this section that h = 1, i.e., the regret function reduces to (1.2). 

3.1. Double exponential jump diffusion. The double exponential jump diffusion process is a Levy process of 
the form 

N t 

(3.1) X t := x + fit + aBt + ^Zi, t>0 

i=i 

where \i G R, a > 0, B = {B t ; t > 0} is a standard Brownian motion, N = {N t ; t > 0} is a Poisson process 
with parameter A > and Z = {Zf, i G N} is a sequence of i.i.d. random variables having a double exponential 
distribution with its density 

(3.2) f(z) := Pn-e v - z l {z<Q} + (1 - p)7] + e-^ z l {z>0} , z G M, 

for some 7y_, r] + > and p G [0, 1]. Here B, N and Z are assumed to be mutually independent. 
The Laplace exponent of this process is given by 

(3.3) V>(/3):=IE M + ^V + Ap^ + fc^-lV /? G M. 

We later see that the Laplace exponent and its inverse function are useful tools in simplifying the problem and 
characterizing the structure of the optimal solution. 

Fix q > 0. There are four roots of ip(P) = q, and in particular we focus on and £2,13 such that 

< < V- < < 00 and ^(-£l,g) = ^(-£2,9) = <?• 
Suppose that the overall drift is denoted by u := E°[Xi], then it becomes 

. \ f P , l ~P 
u = 11 + A h 



n- n+ 



and 



(3-4) 6,0 < " < . > and 6,<? -»• 6,0 as g -)• 

[ > 0, u > j 

for some £x,Q and ^2,0 satisfying 

< 6,0 <V-< 6,0 < 00 and ^(-£i,o) = ^(-6,0) = 0; 
see Figure 1 for an illustration. When u < 0, by (3.4), l'Hopital's rule and £i q ) = q, we have 

0-5) ^^-^ = 4 = ^. 

q tp '(U) u \u\ 

We will see that these roots characterize the optimal strategies; the optimal threshold levels can be expressed in 
terms of and &,q when q > and £2,0 and u when q = 0. 
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(a) U < 



(b) u > 



Figure 1 . Illustration of the Laplace exponent and its roots when the drift u is negative 
and positive. 



Due to the memoryless property of its jump-size distribution, the violation risk and regret can be obtained 
explicitly. The following two lemmas are due to Kou and Wang [20], Theorem 3.1 and its corollary. Here we let 



1,9 



" ^ > and l 2>q :^ " r? " 



6,9 - e 



1,9 



6,9 - 6,9 

where, in particular, when q = and u < 0, 

h,o = "r- > and l 2fi = 6 '° ~ - > 0. 

6,0 6,0 

Notice that + ^2,9 = 1 for every q > 0. 

Lemma 3.1 (violation risk). For every q > a«<i < A < x, we have 



> 0, g > 0, 



Ri q) ( 



-T]-A 



TA) 



(6,9 " 7?-)/i, g e-^- A ) - - ii, q )h, q e~^ x - A) 



In particular, when q = ant/ it < 0, this reduces to 

R W( TA ) = l 2fie -^ A [l - e -^,o(-A) 

Lemma 3.2 (functional associated with the regret when h = 1). For every q > 0, we have 



K 1 



f' T A 

Jo 



- qt dt 



1 



hA2,a ( 1 - e-^~ A A + l 2 M, q ( 1 - e -6.«(-A) 



< A < x. 



Furthermore, it can be extended to the case q = 0, by taking q I via (3.5) and the monotone convergence 
theorem; 

\x - A) + &^!h 



IE [t a ] 



i 

N . 
oo, 



»?-?2,0 



.(1 _ e -6,o(a!-A)\ 



ij u < 0, 
ifu>0. 
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For a general Levy process, Lemma 3.2 can be alternatively achieved by obtaining F x |x e ^ < A j where e q 
is an independent exponential random variable with parameter q > and X t := info< u <tX u is the running 
infimum of X; see Bertoin [5], Kyprianou [22] or Chapter 2 of Surya [39]. In particular, F x |^ eq < ^} admits 
an analytical form when jumps are of phase-type (Asmussen et al. [1]); the above results can be seen as its special 
case. 

Remark 3.1. As in Kou and Wang [20], we assume throughout this section that X contains a diffusion component 
(a > 0). Here we do not consider the case a = because its spectrally negative case is covered in the next section. 
The results for a = can he obtained similarly using the results by Asmussen et al. [1]. 

3.2. Optimal strategy when h = 1. We shall obtain the optimal solution for q > 0. When q = 0, we focus on 
the case when u < because K x 9 = oo otherwise by Lemma 3.2 and the problem becomes trivial as we discussed 
in Section 2. 

Suppose q > 0. By Lemma 3.2, the stopping value (2.3) becomes 



(3.6) 



G(x) = 7E* 



where 
(3.7) 

which satisfy 
(3.8) 



~ qt dt 



Ci„ 



i=l,2 



7 h,q&,q 



and C2 , 



£ C ijq e-^ x , x>0, 
i=l,2 

7 h,q£l,q 



i=l,2 



7 ^, »7- 



and £ Ci i9 



7 



n- 



i=l,2 ■'- * 

The difference between the continuation and stopping values (defined in (2.5)) becomes, by Lemmas 3.1-3.2, 



(3.9) S A (x) 



V- 



(6,, - V -)e-^ (x - A) - v - A - ( 1 - '-^-^ 



+ 



2,9 



-(77- - 6, g )e- ?2 ^ (:r - A) - ,7 - A - (l " e- 6 -«(*- A )) 



< A < x. 



Suppose g = and u < 0. By Lemma 3.2, we have 

G(x 



_7_ 

lul 



x > 0, 



and taking g — > in (3.9) via the monotone convergence theorem (or by Lemmas 3. 1-3.2) 

V- 



(3.10) S A (x) 



(6,0 - V-K^ ~ i^6,o(x - A) 



+ 



2.0 



V- 



-r/_e 



-6,o(x-A)-r;_j4 _ 7 



1 



6,0(1- A) 



< A < x. 



Remark 3.2. We have = for every A > 0, i.e., continuous fit holds whatever the choice of A is. This 

is due to the fact that X has paths of unbounded variation (a > 0). As we see in the next section, continuous fit is 
applied to identify the optimal threshold level for the bounded variation case. □ 
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S'a(A- 



'/- 



We shall obtain the threshold level A* such that the smooth fit condition holds, i.e., 5' A ,(A* +) = if such a 
threshold exists. By (3.9)-(3.10), we have 

\v--tl iq )fo*-V-)e-*- A -*Zi iq Z2, q 

Therefore, on condition that 

' (V- ~ 6,9X6,9 - »7-) > ^6,96,9; 9 > 



g > 

q = and u < 



A > 0. 



(3.11) 



(6,o -»?-)> jiJr6,o, 



g = and u < 

the smooth fit condition 5' A (A+) = is satisfied if and only if A = A* > where 



(3.12) 



A* :-- 



7 £l,g&2,g 



n_ 6 V « 



g>o, 

g = and u < 0. 



If (3.11) does not hold, we have d' A (A+) < for every A > 0; in this case, we set A* = 0. 

We now show that the optimal value function is <p '■= 4>A* (see (2.6)). Suppose q > and A* > 0. Simple 
algebra shows that 

7 1 



6,9 e-^-^)-l -£i, g ( e -6.«<*-A->-i 



(3.13) 5(x):=^(x) 

9 6,9 - 6,g 

This together with (3.6) shows that 

' Ei=i, 2 (^,9 - Ci, g )e-^ x ; 



x > A* 



1, 



where Ci j9 and C2, g are defined in (3.7) and 



(3.14) 



7 6,. 



q 6, 9 - 6,g 

When <7 = and A* > 0, we have 



z^- qA * and L 2 „ := 



x < 0, 



7 6,9 „£ 2 , a A* 



S(x) 



\u\ 



(x - A*) 



I — e -6,o(a!-A*) 
6,0 



x > A*. 



and consequently, 



7 (a* + ^ 2 '°~' ? - h - p-6,o^ _i_ J_fl _ p -6,o(z-^* 
|n| l A + »7-6,o y L e ) + 6,o^ e 



(*) = < 



|u| 
1, 



x > A*, 
< x < A* 
x<0. 



Finally, it is understood for the case A* = (for both q > and q = 0) that 
(3.15) 

where, by Lemma 3.1, 



(x) = ]im(f) E (x) 
s4-0 



Hmi?W(r £ ) = E* [e^ Vs<0, *<oo} 



(?7- ~6,g)(6,g -f\- 
^-(6,9-6,9) 
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This is the expectation of the cost incurred only when it jumps over the level zero. 

When A* > the value function can be attained by ta* , and when A* = it can be approximated arbitrarily 
closely by <f) e (x) (which is attained by r e ) for sufficiently small e > 0. We therefore only need to verify that 
<j>(x) < W [e- qr G(X T )l {T<oo} ] for any r E S. 

We first show that </>(•) is dominated from above by the stopping value G(-). 

Lemma 3.3. We have <p(x) < G(x) for every i€R 

Proof. For every < A < x, 8(/)a(x) / dA = 85a{x) /OA equals 

± [h, g e-^ x - A ^ + hae-^*-*)] [-(v- ~ Ci,M2, q ~ ri-)e- v - A + , q > 

£ [h,0 + l2,oe-^ o{x ~ A) ] [-r?_(6,o - r / „)e""- A + ^6,o] , q = and u < 

Here, in both cases, the term in the first bracket is strictly positive while that in the second bracket is increasing in 
A. Therefore, when A* > 0, A* is the unique value that makes it vanish, and consequently d(j)A(x)/dA > if 
and only if A > A*. On the other hand, if A* = 0, d4>A(x)/dA > for every < A < x. These imply when 
x > A* that by Remark 3.2 

4>{x) < lim 4>a{x) = G(x) and <j)(x) = \im.<j) £ (x) < lim^fx) = G(x) 

when A* > and when A* = 0, respectively. On the other hand, when — oo < x < A*, we have <j){x) = G{x) by 
definition and hence the proof is complete. □ 

Remark 3.3. (1) Suppose q > 0. In view of (3.6), G(-) is bounded from above by j/q uniformly on x E 

(0,oo). 

(2) Suppose u < and fix x > 0. We have G(x) < ¥, x 9 < oo uniformly on q > 0; namely, G(x) is bounded 
by K x 9 < oo uniformly on q E [0, oo). 

(3) Suppose A* > 0. Using the same argument as in the proof of Lemma 3.3, we have <f>(x) < lim e ^o 4>e{x) = 
W [e"' ?6, l{x () <o,9<oo}] < 1 for every x > A*. When < x < A*, using the monotonicity of G(-) and 
continuous fit (see Remark 3.2), we have 4>(x) = G(x) < G(A*) = <t>(A*) < lim^ 4>e(A*) < 1. When 
A* = 0, we have <j>{x) < 1 in view of (3.15). Therefore, <p(-) is uniformly bounded and it only takes 
values on [0, 1]. □ 

Now we show that C<p(x) > q<fi(x) on (0, oo)\{^4*} where C is the infinitesimal generator of X such that 

1 f°° 

(3.16) Cw(x) = fiw'(x) + -a 2 w"(x) + A / [w(x + z) — w(x)] f(z)dz, iGK, 

2 J-oo 

for any C 2 -function w : M — > K. The proof for the following lemma is lengthy and technical, and therefore is 
relegated to the appendix. 

Lemma 3.4. (1) If A* > 0, then we have 

(3.17) £<f>(x)-q<f>(x) = Q, x > A* , 

(3.18) C(f>(x) - q<f>(x) > 0, < x < A*. 
(2) If A* = 0, then (3.17) holds for every x > 0. 
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Lemmas 3.3 and 3.4 show the optimality. The proof is very similar to that of Proposition 4.1 given in the next 
section (see Appendix A.6) and hence we omit it. 

Proposition 3.1. We have 

<t>{x) = inf E x [e- qT G{X T )l {T<oo} ] , x > 0. 

4. Spectrally Negative Case 

In this section, we analyze the case for a general spectrally negative Levy process. We shall obtain the optimal 
strategy and the value function in terms of the scale function for a general h. We assume throughout this section 
that q > 0. The results obtained here can in principle be extended to the case q = on condition that ¥, x 9 < oo 
along the same line as in the discussion in the previous section. The proofs of all lemmas and propositions are 
given in the appendix. 

4.1. Scale functions. Let X be a spectrally negative Levy process with its Laplace exponent 



cf3 + \<J 2 P 2 



(0,oo) 



-/3x 



i + /3xi {0 < x <i } )n(dx), pe 



where c € M, a > and II is a measure on (0, oo) such that Jj Q ^(1 A x 2 )LI(dx) < oo. See Kyprianou [22], 
p. 212. In particular, when 



(4-1) 

we can rewrite 

where 
(4.2) 



/ (1 Ax)Il(dx) < oo, 

V(0,oo) 



-fix 



(0,oo) 



i)n(dx), pe 



H ■= c + 



/ xn(dx). 

7(0,1) 



The process has paths of bounded variation if and only if a = and (4.1) holds. It is also assumed that X is not a 
negative subordinator (decreasing a.s.). Namely, we require fi to be strictly positive if a = 0. 
It is well-known that tp is zero at the origin, convex on M + and has a right-continuous inverse: 

( q := sup{A > : ^(A) =q}, q> 0. 

Associated with every spectrally negative Levy process, there exists a (q-)scale function 

Wis) : R -> E ; g > 0, 

that is continuous and strictly increasing on [0, oo) and satisfies 



1 



If r+ is the first time the process goes above a > x > 0, we have 

W(«)(a;) 



{t+<8,t+<oo} 



WW (a) 



and E a 



'1 



{r+>e,e<oo} 



- Z^(o) 



VK(?)(a) 
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where 

(4.3) Z^(x):=l + q[ W {q \y)dy, x£R. 

Jo 

Here we have 

(4.4) W {q \x) = on (-oo,0) and Z iq \x) = 1 on (-oo,0]. 

We assume that II does not have atoms; this implies that W^) is continuously differentiable on (0, oo). See Chan 
et al. [9] for the smoothness properties of the scale function. 
The scale function increases exponentially; indeed, we have 

C x 

(4.5) W iq) (x) ~ 4ttt\ asx^oo. 

ViCq) 

There exists a (scaled) version of the scale function = {Wq(x); that satisfies, for every fixed q > 0, 

W (q (x) = e~ c " x W {q) (x) , i£i 

and 

Moreover (x) is increasing and as is clear from (4.5) 

(4.6) W <q (x) /* — — asi^oo. 

V (Cq) 

From Lemmas 4.3 and 4.4 of Kyprianou and Surya [24], we also have the following results about the behavior 
in the neighborhood of zero: 

(4.7) 



( „\ , , 0, unbounded variation /„ v , 

W^(0) = ^ ' ^ and (0+) 

\ i , bounded variation J 



oo, cr = and 11(0, oo) = oo 

g+n(o,oo) ^ com pound Poisson 



For a comprehensive account of the scale function, see Bertoin [5, 6], Kyprianou [22] and Kyprianou and Surya 
[24]. See Surya [40] and Egami and Yamazaki [13] for numerical methods for computing the scale function. 

4.2. Rewriting the problem in terms of the scale function for a general h. We now rewrite the problem in 
terms of the scale function. For fixed A > 0, define a random measure 

M^9)(w,B) := [ TAiUl) e- qt l {Xt{w)eB} dt, B(R). 
Jo 



Lemma 4.1. For any u £ £1, we have 
(4.8) / e 



o 



qt h{X t (u))dt = [ M^ q \u,dy)h(y). 
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With this lemma and the property of the random measure, we have 



(4.9) 
where 



e- qt h(X t )dt 



M^\u,dy)Ky) 



^ q \dy)h{y) 



^ x A >i\B):=E x M( A ' q \B) , B G B(R) 

is a version of the q-resolvent kernel that has a density owing to the Radon-Nikodym theorem; see Bertoin [6]. By 
using Theorem 1 of Bertoin [6] (see also Emery [15] and Suprun [38]), we have for every B G B(M) and a > x 

■W( q \x - A)W®(a - y) 



e- qt l {Xt£B} dt 



Bn[A,oo) 



WM(a-A) 



hx>v}W®(x-y) 



dy 



Cq{y _ A) W^ X -A)W (q (a-y) _ (a) , _ } 



'Bn[A,oo) 

Moreover, by taking a f oo via the dominated convergence theorem in view of (4.6), we have 



dy. 



lim 

Bn[A,oo) a ^°° 



_ (q{y _ A) W( q H X -A)W (q (a-y) _ _ 



< q (y-A) w (q)^ x -A)- l {x > y} W( q \x - y) 

'Bn[A,oo) 

where the second equality holds by (4.6). Hence, we have the following result. 



{x>y} ' 
dy 



dy 



Lemma 4.2. Fix q > and < A < x. We have 

( e -< q (y-A) w ( q ) {x _ A) _ i {g£ty}W h\ x - y ))dy, y > A, 

y <A. 



»i A > q Hdy) 



0. 



By (4.9) and Lemma 4.2, we have, for any arbitrary < A < x, 



(4.10) 



TA 



e- qt h(X t )dt 



W (q \x-A) / e-<" y h(y + A)dy- / W^(x - y)h{y)dy, 



and this can be used to express the regret function and the stopping value. This also implies that the integrability 
condition (1.1) is equivalent to f °° e~^ qV h(y)dy < oo. 

Using the g-resolvent kernel, we can also rewrite the violation risk. 



Lemma 4.3. For every < A < x, we have 

rM (j a ) = —W M (x-A) / n(du) ( 1 - e 

Cq J A V 



< q (u-A) 



1 

1 J A 



n(d«) 



A) - Z®(x - «)) . 



By (4.10) and Lemma 4.3, the difference function (2.5) becomes, for all < A < x, 



8 A (x) = R ( i\r A )- 1 E x 



e- qt h{X t )dt 



"1 /* OO 1 /* OO 

(4.11) = -W iq \x -A) H(d«) (l - e -^- A )) - - / n(dn) (z®(x -A)- Z®( 

/*oo f*X 

-7 WW(x-A) / e~^ y h(y + A)dy - / (x - y)h(y)dy 

l Jo Ja 



x — u 
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4.3. Continuous and smooth fit. We now apply continuous and smooth fit to obtain the candidate threshold levels 
for the bounded and unbounded variation cases, respectively. Firstly, the continuous fit condition 6a{A+) = 
requires in view of (4.11) that 



(4.12) 
where 



W<«>(0)*(A) = 



1 (* OO /* OO 

$(A):=— / Il(du) ( 1 - e-^-^) - 7 / e~^ y h{y + A)dy, A>0. 
Gj J A v ' Jo 

Note in this calculation we used the fact that the second term in (4.11) vanishes as x \. A by (4.4). The condition 

(4.12) automatically holds for the unbounded variation case by (4.7), but for the bounded variation case it requires 

(4.13) §(A) = 0. 

For the unbounded variation case, we apply smooth fit. For the violation risk, by using (4.3)-(4.4), we have 



d_ 
dx 



R x q) ( 



TA 



x=A+ 



1 r 00 / 

W (i)\ 0+ ) / n(dn) (1 - e~^ u - 

Cq J A V 



■A) 



For the regret function, by using (4.7) in particular W^(0) = 0, we obtain 



d_ 
dx 



E 1 ' 



TA 



e - qt h{x t )dt 



= lim 

x=A+ X ^ A 



W iqY {x-A) / e~^ y h{y + A)dy- / W iq)> (x - y)h(y)dy 



poo 

= iy (<?) '(0+) / e -t« y h(y + A)dy. 
Jo 

Therefore, for the unbounded variation case, smooth fit requires 

W {q) ' '(0+)*(>l) = 0. 

Consequently, once we get (4.13), continuous fit holds for the bounded variation case and both continuous and 
smooth fit holds for the unbounded variation case (see Figure 4 in Section 5 for an illustration). Because h is 
non-decreasing by assumption, we have 



n(du)e 



-C(u-A) 



7 



/•oo 

/ e-^ y h'(y + A)dy < 0, A > 0. 
Jo 



Hence there exists at most one root that satisfies (4.13). We let A* be the root if it exists and zero otherwise. 
Because linutoo ®{A) < 0, A* = means that $(A) < for all A > 0. 



4.4. Verification of optimality. We now show as in the last section that the optimal value function is (ft := (ft 
(see (2.6)) where in particular the case A* = is defined by (3. 15). When A* > 0, it can be attained by the strategy 
ta* while when A* = it can be approximated arbitrarily closely by r £ with sufficiently small e > 0. Recall that 
(ft{x) = 5(x) + G(x) (with 5 := 5a*) for x > A* and that G{x) in (2.3) for x > can be expressed in terms of the 
scale function by taking the limit A J, in (4.10). The corresponding candidate value function becomes both for 
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A* = and for A* > 
(4.14) 

Six 



0, 

+ 7 



-w^{x - A*) / n(d«) (i - e -<^- A *)) - - / n(dtt) (z {q \x - A*) - Z^\x - u) 

,q J A* V ' Q J A* V 



W {q \x) / e-^ y h{y)dy-W ( - q) (x- A*) / e~^ y h(y + A*)dy - / - y)h(y)dy 

Jo Jo Jo 



for every x > 0. By definition, <^(x) = 1 for every x < 0. In particular, when A* > 0, we can simplify by using 
(4.13), 



<f,{x) = - J y°° n(dn) (z®(x - A*) - Z^(x - «)) 



+ 7 



/■oo /"A 

iy (9) (x) / e-^ y h{y)dy- / W (g) (x - y)h(y)dy 
Jo Jo 



x>0. 



These expressions for the candidate value function 6 are valid not only on (A*, oo) but also on (0, A*] thanks to 
(4.4) and continuous fit. 

In order to verify that 6 is indeed optimal, we only need to show that (1) 6 is dominated by G and (2) C6(x) > 
q&{x) on (0, oo)\{^4*} where 

Cf(x) = cf{x) + -a 2 f"(x) + J [f{x -z)- f(x) + f'(x)zl {0<z<1} ] U(dz), 

poo 

Cf(x) = pf'{x) + / [f(x -z)- f(x)) U(dz), 
Jo 

for the unbounded and bounded variation cases, respectively; see (4.2) for the definition of fi. The former is proved 
in the following lemma: 

Lemma 4.4. We have 6(x) < G(x) for every iei 

Now we recall that the processes | e -«(* AeAr « ) W^ q \X tAdAT +); t > o} and [e~^ tA9AT ^Z^(X tA0AT +); t > o} 
are -martingales for any < x < a; see page 229 in Kyprianou [22]. Therefore, we have 

(4.15) CW ( - q \x) = qW^(x) and £Z®(x) = qZ®(x), x > 0. 

We take advantage of (4.15) to show the following lemma. 

Lemma 4.5. (1) We have 

C6[x) - q6(x) = 0, x>A*. 

(2) If A* > 0, we have 

C6(x) - q6{x) > 0, < x < A*. 
Finally, Lemmas 4.4-4.5 are used to show the optimality of 6. 
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Proposition 4.1. We have 

4>{x) = iuf E* [e-« r G(JT r )l {r<oo} ] , x > 0. 
5. Numerical Results 

We conclude this paper by providing numerical results on the models studied in Sections 3-4. We obtain optimal 
threshold levels A* for (1) the double exponential jump diffusion case with h = 1, and for (2) the spectrally 
negative case with h in the form of the exponential utility function. We study how the solution depends on the 
process X. We then verify continuous and smooth fit conditions for the bounded and unbounded variation cases, 
respectively. 

5.1. The double exponential jump diffusion case with h = 1. We evaluate the results obtained in Section 3 
focusing on the case h = 1. Here we plot the optimal threshold level A* defined in (3.12) as a function of 7. The 
values of and £2,9 are obtained via the bisection method with error bound 10~ 4 . 

Figure 2 shows how the optimal threshold level changes with respect to each parameter when q = 0.05. The 
results obtained in (i)-(iv) and (vi) are consistent with our intuition because these parameters determine the overall 
drift u, and A* is expected to decrease in u. We show in (v) how it changes with respect to the diffusion coefficient 
cr; although it does not play a part in determining u, we see that A* is in fact decreasing in a. This is related to the 
fact that, as a increases, the probability of jumping over the level zero decreases. 

5.2. The spectrally negative Levy case with a general h. We now consider the spectrally negative case and 
verify the results obtained in Section 4. For the function h, we use the exponential utility function 

h{x) = 1 - e~ px , x>0. 

Here p > is called the coefficient of absolute risk aversion. It is well-known that p = —h"(x)/h(x) for every 
x > 0, and, in particular, h = 1 when p = 00. We consider the tempered stable ( CGMY) process and the variance 
gamma process with only downward jumps. The former has a Laplace exponent 

$(P) = c(3 + CX a T(-a) { f 1 + £) " - 1 - ^ J 

and a Levy density given by 

n ( d2; ) = C ^l+^ 1 {x>0} dx 

for some A > and a < 2; see Surya [39] for the calculations. It has paths of bounded variation if and only if 
a < 1. When a = 0, it reduces to the variance gamma process; see Proposition 5.7.1 of Surya [39] for the form 
of its Laplace exponent. We consider the case when a = 0. The optimal threshold levels are computed by the 
bisection method using (4.13) with error bound 10~ 6 . 

Figure 3 shows the optimal threshold level A* as a function of 7 with various values of p. We see that it is indeed 
monotonically decreasing in p. This can be also analytically verified in view of the definition of &(A) because h is 
decreasing in p for every fixed A, and consequently the root A* must be decreasing in p. This is also clear because 
the regret function monotonically decreases in p. 
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Figure 2. The optimal threshold level A* with respect to various parameters: the pa- 
rameters are /i = —1, a — 1, rj- — 1.0, r] + = 2.0, p = 0.5, A = 1 unless they are 
specified. 
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p=0.01 

- p=0.02 
p=0.05 
p=0.1 




(a) tempered stable (unbounded variation) 




(b) tempered stable (bounded variation) 




(c) variance gamma 

Figure 3. The optimal threshold level A* with various values of coefficients of absolute 
risk aversion p: (a) tempered stable (unbounded variation) with a = 1.5, A = 2, C = 0.05 
and c = 0.05, (b) tempered stable (bounded variation) with a = 0.8, A = 2, C = 0.075 
and c = 0.05, and (c) variance gamma with with A = 2, C = 0.075 and c = 0.05. 



5.3. Continuous and smooth fit. We conclude this section by numerically verifying the continuous and smooth 
fit conditions. Unlike the optimal threshold level A*, the computation of the value function involves that of 
the scale function. Here we consider the spectrally negative Levy process with exponential jumps in the form 
(3. 1) with p = 1 and a > 0. We consider the bounded variation case (a = 0) and the unbounded variation case 
(<r > 0). We also set h = 1. This is a special case of the spectrally negative Levy process with phase-type jumps, 
and its scale function can be obtained analytically as in Egami and Yamazaki [13]. In general, scale functions 
can be approximated by Laplace inversion algorithm by Surya [40] or the phase-type fitting approach by Egami 
and Yamazaki [13]. One drawback of the approximation methods of the scale function is that the error tends to 
explode as x gets large (see (4.5)). Because our objective here is to accurately verify the continuous and smooth 
fit conditions, we use an example where an explicit form is known. Notice that the threshold level A* can be 
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computed independently of the scale function, and hence one can alternatively approximate the value function <fr 
by simulation. 

Figure 4 draws the stopping value G as well as the value function for both the bounded and unbounded 
variation cases. The value function is indeed continuous at A* for the bounded variation case and differentiable 
for the unbounded variation case. It can be seen that G indeed dominates <fi. While G is monotonically increasing 
on (0, oo), <j) is decreasing for large x and is expected to converge to zero as x — > oo. 




(a) continuous fit 



(b) smooth fit 



Figure 4. Illustration of <p (solid) and G (dotted) for (a) the bounded variation case with 
(j = and fj, = 0.3 and (b) the unbounded variation case with a = 0.5 and [i = 0.175. 
Other parameters are q = 0.05, h = 1, X = 0.5, 77 = 1 and 7 = 0.04. 



Appendix A. Proofs 



A.l. Proof of Lemma 3.4. We shall prove for the case q > and then extend it to the case q = and u < 0. We 
first prove the following for the proof of Lemma 3.4. 

Lemma A.l. Fix q > and x£l. Suppose that a function w : M — ?■ M in a neighborhood of x > is given by 

){x) = k + ^2 k ie~^' qX 



w( 



i=l,2 



for some k, k± and ki in K. Then we have 
q(w(x) — k) = Cw{x) — A 



+ z)f(z)dz- \ k+J2 k * 
i=l,2 



m- , 0--p)v+ 



l.q 



Proof. Because tp(— £i )(J ) = tp(— £2,5) = q, we have 

q(w(x) - k) = ^2 ki^{-£i, q )e~ il 



i=l,2 
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Moreover, the right-hand side equals by (3.3) 



i=l,2 



(1 - g>7+ _ 1 



*7- - 6,? + 6 



-A(w(x) - k) + V" fe 

i=l,2 

-a 2 w"(x) + fiw'(x) — Xw(x) + Xk + A fe 



(1 ~ g)ffj 
V- ~ 6,9 ' V+ + 6,9 



+ 



i=l,2 



(1 ~ P)^ 
»7- - 6,9 ' ??+ + 



+ 



Cw(x) — A 



/ u; (a + z)/(z)dz - [ fc + ^2 k 

J — OO \ • i n 



8=1,2 



+ 



»7- - 6,9 r}+ + 6,9 



as desired. 



□ 



Proof of Lemma 3.4 when q > 0. (i) Suppose A* > 0. By Lemma A. 1 above, we have C<f>{x) — q<ft{x) equals, for 
every x > A*, 



(A.1) 



A 



/oo 
<f>(x + z)f(z)dz- J2( L i,q-°i,q) 

-°° i=l,2 



+ ~ ) e -e 



v- - 6,g »?+ + 6 



t,9 



and, for every < x < A*, 



(A.2) 



-(g + A^ + A 



/oo 
0(x + z)f(z)dz + ^ Ci,, 
■°° i=l,2 



(1 ~ P)^ 

»7- - 6,g ' »7+ + 6,9 



+ 



by using (3.6) and (3.14). 

Proof of (3.17). We only need to show (A.l) equals zero. Notice that the integral can be split into four parts, 
and, by using (3.7)-(3.8), we have 



poo 

/ <f>(x + z)f(z)dz = y2 ( L *,q ~ C i,q 

J ° i=l,2 



0--p)v+ ( 
v+ + 6,9 ' 







(x-A*) 
(x-A*) 



+ Z)f(z)dz = (Li,q ~ C iy q)— ,--'■<' 



i=l,2 



V~ ~ 6,9 



Pit c -v-(x-A*)-ti.aA* 



i=l,2 



v- - 6,9 



(x + z)f(z)dz 



2=1,2 



i,9 



PV- c - v _(x-A*)-Z t . A* 



v- - 6,9 



(/)(x + z)f(z)dz = pe~ ri - x . 
Putting altogether, (A.l) equals 



J i,9 



2=1,2 



'■9 



+ — e 
Q 



+ pe 



pe 



-Tj—X 



7 



6,96,. 



(V- ~ 6,9) (6,9 -v-) 
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and this vanishes because of the way A* is chosen in (3. 12). 

Proof of (3.18). We shall show that (A.2) is decreasing in x. Note that J™ <p(x + z)f(z)dz can be split into 
four parts with 



OG 

I 

A*-x 
A*—x 



+ z)f(z)dz = (1 -p) (Li, q ~ C i>g )-^—e-^ A *-^<> A * , 

69 + V+ 



(x + z)f(z)dz = (1 — p) 



i=l,2 

7 



o 



1-e 



- v+ (A*-x) 



i=l,2 



r l+ I e -a^i, s _ e -r?+(A*-x)-^,,A* 



9 6,9 + ^7+ 



+ z)/(z)d2; 



p E ^ 

i=l,2 



0(a; + z)/(z)dz = ' ? - :r , 
and after some algebra, we have 

/oo 
<j>{x + z)f(z)dz + Ci, q 
-oo ■ i o 



i=l,2 



PV- , (1 -p)ri+ \ 



J,9 



+ pe _,? - x + (1 -p)e 



ry+(x-A*) | _7 



+ E^ 

£=1,2 



6,g + »7+ 



which by (3.14) equals to 



7 



7 6,96,9 / 1 



q 6,9-6,9 9 V6,9 + r ?+ 6,9+^4 



1 



Hence we see that (A.2) or C4>(x) — q4>(x) equals 



-7 + A 



pe v ~ x - (1 - p)e 



-ri+(A*-x) 6,g6,9 



7 



<? (6,9 + r /+)(6,9 + ^?+). 

and is therefore decreasing in x on (0, ^4*). We now only need to show that lim^. (£4>(x) — q4>(x)) > 0. 
For every x > A* , 

7 6,96,9 



S"(x) 

and hence, after taking x I A*, 



q 6,9 - 6,g 



6,9^ 1 ^ A * ) -6,ge- 6 - (x - A * ) 



6"(A*+) = -^6,96,9 < 0. 

Consequently, by the continuous and smooth fit conditions 5(A*+) = 5'{A*+) = and the definition of C in 
(3.16), we have 

lim (£6(x) — q4>{x)) > lim (C(j)(x) — qMx)) = 0, 
x^A* xlA* 



as desired. 
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(ii) Suppose A* = 0. By (3.15), we have 



C (e-^-" x - e -^ x ) , x > 



where 



c . = (V- ~ 6,9X6,9 -v-) 
r?-(6,g-6,g) 

By Lemma A.l above, we have C<f)(x) — q4>(x) equals, for every x > 0, 



/oo 
<P(x + z)f(z)dz- C 
■°° i=l,2 



+ 



?7- - 6,9 ??+ + 6. 



and this vanishes after some algebra. 



□ 



Proof of Lemma 3.4 when q = 0. We extend the results above to the case q = and u < 0. Solely in this proof, let 
us emphasize the dependence on q and use A* 4> q (-), G q {-), S q (-), and 5a,q{-) with a specified discount rate q > 0. 
We shall show that C4>q{x) = lim g _»o C(j> q {x) for all x > 0. First notice that ^4* — > Aq as g — > 0. 
Clearly, G q (x) — > Gq(x) by the monotone convergence theorem. Furthermore, we have 



G' q (x) = -^6, g 6, 9 £ /;,,< - ^> ^ 
J gr/_ ' |u| 



1 + ^°-^- e -e 
V- 



2,QX 



G'o(x), 



940. 7 



(6,o - V-A&fl _ €: 



2,0a; 



Gg(x) 



Fix x > Aq. We suppose Aq = and focus on g G [0, go] with qo > sufficiently small such that ^4* = for 
all < q < qo. We have, by applying the monotone convergence theorem on (3.15), <p g (x) — > 4>o(x) as q — > 0. 
Furthermore, 



^2 -ii,qX, 



y r?-(6,g-£i,<z) V 



9-1-0, 



-6,0a; 



94.0, 



?2,0 e 



Sg(x). 



Suppose Aq > and focus on g G [0, go] with go > sufficiently small such that x > A* for all < g < go. Note 
that 

\5 g (x) - 6 (x)\ < \5 q (x) - 6a*,o(x)\ + \Sa*,q( x ) ~ s o( x )\ 

where on the right-hand side the former vanishes as q — > by the monotone convergence theorem in view of (2.5) 
and the latter vanishes because A* — > Aq; hence 5 q (x) — > 5o(x) as q — > 0. Furthermore, 

7 6,96,9 



Sg(x) 



q 6,g - 6,9 - 

7 6,96,9 

q 6,9 - 6,9 



_ e -6, 9 (2-A*) +e -&, q (x-A*) 



94Q ; 7 

lul L 



1 _ e -6,o(a:-^o) 



5q{x), 
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In summary, we have 4> q (x) — > 4>o{x), 4>' q ( x ) ~^ 4>b( x ) an d 4>q( x ) ~^ ftoix) as Q — > f° r every x > 0. 
Moreover, by Remark 3.3-(3), via the dominated convergence theorem, 

/oo poo 
(f>o(x + z)f(z)dz = lim / 4> q (x + z)f(z)dz, x>0. 
-oo 1^°°J-oo 

Consequently, lim 9 _>o {£4>q{ x ) ~ Q4>q{ x )) = C4>o( x )- This together with the result for q > shows the claim. □ 

A.2. Proof of Lemma 4.1. Because h is continuous on (0, oo), it is Borel measurable. Hence there exists a 
converging sequence of simple functions (/i^) n eN increasing to h in the form 



l(n) 



h^(y):=J2t\veB n , i} , n>l, 



i=i 



for some Z : N -> N, {6^ n) ;n > 1,0 < % < l(n)} and Borel measurable sets {B n ^n > 1,0 < i < l(n)}. See 
page 99 of Cinlar and Vanderbei [11]. 

Then the right-hand side of (4.8) is, by the monotone convergence theorem, 

r r Kn) 

/ M^\co,dy)h(y) = / M^(u,dy) lim J2t\yeB ni} 

= lim ^bf / M( A >«\u,,dy)l {yeBni} = lim ^b^ (u, B n ,). 

1=1 JK i=l 



This is indeed equal to the left-hand side of (4.8) because, by the monotone convergence theorem, 



-*ft(X t (a;))dt = / e"^ lim ^6^1 {XtMeBni} dt 
1/0 i=l 

i=l J0 i=l 



as desired. 



A.3. Proof of Lemma 4.3. Let iV(-, •) be the Poisson random measure for —X and X_ t :- 
t > 0. By the compensation formula (see, e.g., Theorem 4.4 in Kyprianou [22]), we have 



min < s <t X s for all 



R^(t A )=E* 
= E a 



/■CXJ f'OO 

/ / iV(dt,dn) e -"*l {Xt __ u < .x t >A} 
Jo Jo 

/ e~ qt dt / n(du)l {Xt __ u < 0i x t >A} 
Vo Jo 

/ E(d«) / dt [e- qt F x {X t - < u,X t _ > A}] 
Jo Jo 



n(d«) / dt [ e - g *p x {Xi_ < li, T A > t}} 

Jo 
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By using the (/-resolvent kernel that appeared in Lemma 4.2, we have for u > A 

' e -C*(v- A )W®(x -A)- W {q) (x - y) 
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dt [e-^iXt- <u,t a > t}] = / dy 

ii J A 

u—A 



dz 







<" z W {q) (x -A)- W iq) (x -z-A) 



c 

and it is zero on < u < A. Substituting this, we have 

1 



u—A 



^W^(x — A) (l — e- c " (u - A) ) - I dzW {q) (x -z-A), 
^ ' Jo 



Ri q) (r A ) 



n(d«) |-iy (9) (^ - A) (l - e-^ u ~ A ^ - J d-AV" l! (. r - A 



u-A 



By (4.4), we have 

Rti)( TA ) = ^W^(x - A) r n(dn) fl - e -^ u ~ A A 

Cq J A ^ ' 

f*x pu— A p oo /*x — A 

- \ n(du) / dzW iq \x - A - z) - / n(du) / dzW {q) (x- A- z) 

J A JO Jx Jo 

= -W {q) (x-A) / Il(du) ( 1 - e~^ {u - A) ) 

Cq J A ^ ' 

--/ n(du) [z {q \x-A) -Z^( x — u)\ / H(du) (z^(x — A) — 

= ±w® (x-A) J U(du) (l - e -^ u ~ A ^ - ± J n(d«) (Z® (x-A)- (x - u)) 



as desired. 



A.4. Proof of Lemma 4.4. Fix < A < x. We have 



d_ 
dA 



R^(t a ) = (w {q) (x -A)- ^W^'tx - A)j J°°n(du) (l - e -^ u ~ A ^ 

1 poo . s 

= W c (x - A)e^ x - A 1 / n(dn) (1 - e -^"- A ) 

Cq J A V 



and because we can write 

<"TA 



we have 



d_ 
dA 



e- qt h(X t )dt 



o 



e (qX W (q (x - A) / e-^ y h(y)dy- / W {q) (x - y)h(y)dy, 
J A J A 



TA 



e~ qt h(X t )dt 



(x - A)e^ x / e-^ y h(y)dy - W Cq (x - A)e^ {x - A) h(A) + W {q) (x - A)h(A) 
J A 



-W; (x - A)e^ {x - A ^ / e~^ y h(y + A)dy. 
Jo 
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Summing up these, we obtain 

^ A (x) = -W' Cq {x - A)e^ x - A ^(A). 
Here, because (x — A) > and $>(A) is decreasing in A and attains zero at A*, we have 

(A.3) ^6 A (x)>0^$(A)<0^A>A*. 
Now suppose A* > 0, we have 

4>(x) = G(x)+5 A *{x) < G(x) + limS A (x) < G(x), x > A* 

A^x 

where the last inequality holds because continuous fit holds everywhere for the unbounded variation case and 
because lim^^. 5 A (x) = W^ q ' (0)<&(x) < by (A.3) for the bounded variation case (by noting that x > A*). The 
case with A* = holds in the same way by the definition that <fi(x) = lim £ ^ 4>s(%)- Finally, because = G on 
(—oo,A*], the proof is complete. 

A.5. Proof of Lemma 4.5. (1) When x > A*, <j) is defined in (4.14). Let 4> be defined such that 4>{x) = 4>{x) for 
all x > and 6(x) = for all x < 0. We obtain 



1 / roo poo \ 

- / U(du)CZ^( x — u) — q / n(du)Z^(a -u))= -IL(x, oo). 



£<fi(x) — qMx) 

q 

Here the first equality holds by (4.15) and because the operator C can go into the integrals thanks to the fact that 
Z( q ' is C 1 everywhere and C 2 on M\{0} for the unbounded variation case and it is C° everywhere and C 1 on 
M\{0} for the bounded variation case, and also to the fact that x > A*. The second equality holds by (4.4). 
Because <p{x) = <p{x) and C(f){x) — C<j)(x) = U(x, oo) for every x > 0, we have the claim. 
(2) Suppose < x < A*. We can write 

poo px 

(A.4) cP(x) = -/W iq) (x) e-^ y h(y)dy--f W {q) (x - y)h(y)dy + L(x), x £ (-oo, A*)\{0} 

Jo Jo 

where L(x) = lr x < \ for every i£i After applying (C — q), the first term vanishes thanks to (4.15). For the 
second term, by integration by parts, 



q 



W {q \x - y)h{y)dy = \h(y)(Z^ (x) - Z^ (x - y))} V * - I ti(y)(Z^(x) - Z^(x - y))dy 



o 



y=o Jo 

h(x){Z^(x) - 1) - / h'(y)(Z (q \x) - Z^\x - y))dy 



X 







h(x){Z^(x)-l)-(h(x)-h(0))Z^{x) + / h'(y)Z^(x-y)dy 

Jo 

-h(x) + h(Q)Z {q \x) + I h'(y)Z ( - q \x-y)dy 

Jo 

pM 

-h(M) + h(0)Z^(x)+ / h'{y)Z^ q \x-y)dy, 



o 
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for any M > x where the last equality holds because Z^ q \x) = 1 on (—00, 0]. The operator (£ — q) can again go 
into the integral thanks to the smoothness of Z^ q > as discussed in (1) and we obtain 

(£ - q) (£ W®(x - y)h(y)dy\ = i(£ - q) f—h(M) + J h'(y)Z^(x - y)dy 



= h(M) + - / h'(y)(C - q)Z^(x - y)dy = h(M) - / h'(y)dy = h(x), 
Q Jo Jx 

where the second to last equality holds by (4.15). For the last term of (A.4), we have 

/•oo 

{C-q)L(x)= / (L(x-u)-L(x))U(du) = U(x,oo). 
Jo 

Putting altogether, we have noting that x < A* and h is increasing, 



(C — q)4>{x) = n(x, 00) — jh(x) 
> J U(du) (l-e 



-t q (u-A*) 



yh(x) 



/oo poo 
n(dn) (l - e -^(«- A *)J - 7 ( q J e -^ y h{x)dy 

n(dn) (l - e -^ u ~ A *^ - -/Cq J e~^ y h(A* + y)dy, 



which is zero because A* satisfies (4.13). This completes the proof. 



A.6. Proof of Proposition 4.1. Due to the discontinuity of the value function at zero, we need to proceed carefully. 
By (2.4) and Lemma 4.4, we must have 1 = 0(0) > 0(0+). 

We first construct a sequence of functions n (-) such that (1) it is C 2 (resp. C 1 ) everywhere except at A* when 
a > (a = 0), (2) 0„(x) = <j>{x) on x G (0, 00) and (3) n (x) f 4>( x ) pointwise for every fixed x G (— 00, 0) 
(with linin^oo n (O) = 0(0+) < 0(0)). It can be shown along the same line as Remark 3.3-(3) that 0(-) is 
uniformly bounded. Hence, we can choose so that n is also uniformly bounded for every n > 1. 

Because 0'(x) = 4>' n (x) and0"(x) = 4>n( x ) on (0, oo)\{A*}, we have (£ — q) (0 n — 0) (x) < for every fixed 
x G (0, oo)\{t4*}. Furthermore, by Lemma 4.5 



f e-i s ((£ - q)<t> n (Xs-))ds\ >E X \ T e~" s ((£ - g)(0„ - 0)(X s _))d. 
Jo i Uo 



(A.5) E x 
Here, the last lower bound is obtained because 



> —00, r G S. 



W 



e-^((£-g)(0 n -0)(AV))ds 



> -KE X 



[ e - qs U(X s ^,oo)ds 
Jo 
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where K < oo is the maximum difference between <j) and <p n . Using N as the Poisson random measure for — X 
and X as the running infimum of X as in the proof of Lemma 4.3, we have by the compensation formula 



E :t ' 



E* 

E x 
E x 

E x 



rU poo roo 

/ e- qs U{X s .,oo)ds =E X 
Jo J Uo Jo 



h 



'0 Jo 

roc poo 



{9> s ,«>x 3 _}n(du)ds 

f'OG f'OO 

/ / e ' qsi {x s >o,«>x s _ } n(du)d s 

Jo Jo 



JO 



e i -{X g <0,9<oo} 



~ qs l{x s _>o,u>x^}N(ds,du) 

< 1. 



By (A.5), we have uniformly in n 
(A.6) 



-qs 



\(C-q)(4> n -4>)(X s .)\ds 



< oo, 



-qs 



\{£-q)((f> n -(f>){X s _)\ds<<x, 



a.s. 



We remark here for the proof of Proposition 3.1 that, in the double exponential case with q = 0, there also exists a 
finite bound because the Levy measure is a finite measure and < oo by assumption. 

Notice that, although <p n is not C 2 (resp. C 1 ) at A* for the case a > (the case of bounded variation), the 
Lebesgue measure of the set where (j) n at which X = A* is zero and hence tp'^A*) (4>' n (A*)) can be chosen 
arbitrarily; see also Theorem 2.1 of [30]. By applying Ito's formula to je -9 ^ 7 ^ <p n (X t /\e)', t > 0}, we see that 



(A.7) 



C ptAO -\ 

le-^ tA9 Un(X tAe ) - e-« s ((C-q)MXs-))ds; t > o| 



is a local martingale. Suppose {o^; k > 1} is the corresponding localizing sequence, namely, 

rtA9/\o k 



K 1 



-q(tA0Acr h ). 



n{Xt/\0/\a kJ 



l {x) + E x 



-qs 



({£-q)</> n (Xs-))ds 



k > 1. 



Now by applying the dominated convergence theorem on the left-hand side and Fatou's lemma on the right-hand 
side via (A.5), we obtain 



- q{tA9 Un(X tAe )] >^ n (x)+E 2 



tA6 



-qs 



((C-q)MX s -))ds 



Hence (A.7) is in fact a submartingale. 

Now fix r € S. By the optional sampling theorem, for any M > 0, 



E J 



- 9 (rAM), 



AM 



>0 n (x)+E 3 



tAM 



-qs 



((C-q)MX s -))ds 



-AM 



-qs 



({£-q)0(X s _))ds 



+ E X 



-AM 



-qs 



((C-q)((t> n -4>)(X^))ds 



where the last equality holds because the expectation can be split by (A.5). Applying the dominated convergence 
theorem on the left-hand side and the monotone convergence theorem on the right-hand side (here the integrands 
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in the two expectations are positive and negative, respectively) along with Lemma 4.5, we obtain 



(A.8) 



E x [e 



t (x r )i {T<oo} ] > Mx) + ^ x 



-qs 



((£-g)(0„-0)(X s _))d S 



For the left-hand side, the dominated convergence theorem implies 



(A.9) lim E x [e-^ n (X T )l {T<oo} ] 



E x 



e qT lim n (X r )l {r<oo} 



= E* [e-i T (<j>(X T )l {XTm + 0(O+)l {XT=o} )l {r <oo}] < [e- qT ^{X r )l {T<oo} ] , 

which holds by equality when a = because X creeps the level zero only when a > 0. For the right-hand side, 
by (A. 6), 



(A. 10) lim E x 

n— >oo 



-qs 



{(C-q){4> n -4>){X s -))ds 



= E X 


[L 







e~ qs hm((C-q)(cP n -cP)(X s _))ds 



Here, for every P x '-a.e. we!!, because X s -(lj) > for Lebesque-a.e. s on (0, t(uj)) 

-<t>)(X s „(u)-u) = 0. 



lim ((£ - g)(^ n - (w)) = / n(d«) lim (</ 

Hence (A. 10) vanishes. 

Therefore, by taking n — )• cxd on both sides of (A.8) (note 4>(x) = <j>n( x ) f° r an Y x > 0)> we nave 

<f>(x) < E x [e^ T ^(X T )l {T<oo} ] < E x [e-« T G(X T )l {T<OQ} ] , r G 5, 

where the last inequality follows from Lemma 4.4. Finally, the stopping time ta* attains the value function <p when 
A* > while t £ and </> e approximate by taking e sufficiently small when A* = 0. This completes the proof. 
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